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ABSTRACT 

We analyze the angular distribution and the orbital rotation directions of a sample of carefully-selected satellite 
galaxies extracted from the Sloan Digital Sky Survey (SDSS). We also study these statistics in an A/-body sim- 
ulation of cosmological structure formation set within the ACDM paradigm under various assumptions for the 
orientations of disk angular momenta. Under the assumption that the angular momenta of the disks are aligned 
with the angular momenta of the inner regions of their host dark matter halos, we find that the fraction of sim- 
ulated satellite halos that exhibit prograde motion is / prog w 0.55-0.60, with larger satellites more likely to be 
prograde. In our observational sample, approximately 60% of the satellites exhibit prograde motion, a result that 
is broadly consistent with the simulated sample. Contrary to several recent studies, our observational sample of 
satellite galaxies show no evidence for being anisotropically distributed about their primary disks. Again, this 
result is broadly consistent with our simulated sample of satellites under the assumption that disk and halo angular 
momenta are aligned. However, the small size of our observational sample does not yet allow us to distinguish be- 
tween various assumptions regarding the orientations of disks in their halos. Finally, we assessed the importance 
of contamination by interlopers on the measured prograde and retrograde statistics. 

Subject headings: galaxies:satellites — galaxies:angular momentum — galaxies:formation — galaxies:evolution 



1. INTRODUCTION 

In the standard hierarchical, cold dark matter (CDM) sce- 
nario of cosmological structure and galaxy formation (e.g. 
White & Rees 1978; Blumenthal et al. 1984), large systems are 
generally formed through the continual merging of compara- 
bly early-forming, low-mass objects in the hierarchy. In this 
context, the halos of large galaxies are formed by smaller ha- 
los that are accreted by the larger host and disrupted by tidal 
interactions. As the massive galactic halo is assembled, other 
halos and their associated galaxies may become gravitationally 
bound to the host and orbit in the host potential before possibly 
being incorporated into the central host galaxy or being strongly 
affected by tides (e.g. Taffoni et al. 2003; Hayashi et al. 2003; 
Zentner & Bullock 2003; Kravtsov et al. 2004b; Kazantzidis 
et al. 2004b; Taylor & Babul 2004; Zentner et al. 2005a). 

The distribution of satellite galaxies about their primary 
galaxies has received significant attention, both theoretically 
and observational. A number of studies have focused on the 
radial distributions of satellite halos and satellite galaxies (e.g. 
Taylor et al. 2003; Zentner & Bullock 2003; Kravtsov et al. 
2004b; Diemand et al. 2004; Gao et al. 2004; van den Bosch 
et al. 2005; Zentner et al. 2005a; Nagai & Kravtsov 2005); how- 
ever, the angular distribution of satellites has become a subject 
of a great deal of recent activity. In an early study, Holmberg 
(Holmberg 1969), found that satellites of spiral galaxies with 
projected separations r p < 50kpc are preferentially located near 
the short axes of the projected light distributions of their host 
galaxies. In other words, they are preferentially located near 
the poles of the primary disk galaxies. Zaritsky et al. (1997a) 
found a statistically-significant anisotropy, similar to that ad- 
vocated by Holmberg (1969), at larger projected separations 
(200kpc< r p < 500kpc) and, in a more recent study, Sales & 
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Lambas (2004) found evidence for the preferential alignment 
of satellites along the minor axes of their primaries in the Two 
Degree Field Galaxy Redshift Survey (e.g. Colless et al. 2001). 
However, the observational status of the so-called Holmberg Ef- 
fect is unclear. Brainerd (2004) studied satellites in the Sloan 
Digital Sky Survey (SDSS, e.g. York et al. 2000; Strauss et al. 
2002) and found evidence for the opposite correlation of satel- 
lite position with the major axis of the light distribution of the 
primary galaxy. 

Kroupa et al. (2005) recently reiterated that the satellite 
galaxies of the Milky Way (MW) are distributed anisotropically 
(see also Lynden-Bell 1982; Majewski 1994; Hartwick 1996, 
2000; Mateo 1998; Grebel et al. 1999; Willman et al. 2004) 
and argued that this observation presents a challenge for the 
standard CDM paradigm. Zentner et al. (2005b, see also Kang 
et al. 2005; Libeskind et al. 2005; Agustsson 2005;) showed 
that MW-sized dark matter halos contain satellite halos or sub- 
halos that are distributed anisotropically about their central host 
halos. Subhalos are located preferentially near the long axes of 
their host dark matter halos, so that the anisotropy of MW satel- 
lites can be understood if the rotation axis of the MW disk is 
closely aligned with the long axis of the triaxial host dark mat- 
ter halo of the MW. In order to test this alignment conjecture, it 
is necessary to examine a large sample of satellites about disk 
galaxy primaries. 

Another interesting way to relate the properties of satellite 
galaxies to disk primaries is through the sense of rotation of the 
satellites relative to the disks, if any. On the theoretical side, the 
prediction for the number of prograde and retrograde satellites 
is unclear and there are many aspects to consider. In the sim- 
plest picture of disk galaxy formation, baryons in halos start 
off sharing the angular momentum distribution of their halos 
and, on average, conserve angular momentum as they cool and 
condense (e.g. Fall & Efstathiou 1980). This implies that the 
poles of disk galaxies should be collinear with the net angular 
momentum vectors of their host halos which, in turn, are gener- 
ally aligned with the halo minor axes (e.g. Warren et al. 1992; 
Porciani et al. 2002; Faltenbacher et al. 2005). Vitvitska et al. 
(2002, see also Mailer et al. 2002) proposed that the angular 
momenta of dark matter halos are built through mergers with 
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infalling satellite halos. If the accretion of satellites is a ran- 
dom process and the acquisition of angular momentum is a ran- 
dom walk as in the Vitvitska et al. (2002) model, this scenario 
seems to favor neither net prograde nor net retrograde satel- 
lite rotation because the angular momentum of the halo at disk 
formation may bear little resemblance to the angular momenta 
carried by the late-accreting satellites. However, the accretion 
of satellite halos generally occurs along a preferred direction, 
so that the accretion directions of satellites are strongly corre- 
lated (e.g. Knebe et al. 2004; Zentner et al. 2005b), the detailed 
consequences of which are unclear. 

There are also dynamical processes that may affect the pro- 
grade and retrograde fractions of satellites subsequent to satel- 
lite accretion. Quinn & Goodman (1986) explored the idea that 
enhanced dynamical friction for satellites on nearly coplanar, 
prograde orbits about the disk could lead to enhanced destruc- 
tion of prograde satellites but found this process to be fairly 
inefficient. Penarrubia et al. (2002) extended this to include 
the effect of an oblate halo about the disk; however, their re- 
sult was that preferred destruction is inefficient for satellites 
further than ~ 50kpc from the primary because the orbital de- 
cay times are large at these radii (see also Zaritsky & Gonzalez 
1999). Moreover, CDM halos generally tend to be more pro- 
late rather than oblate (e.g. Jing & Suto 2002; Bullock 2002; 
Kazantzidis et al. 2004a). Another dynamical argument relates 
to the heating of galactic disks. Velazquez & White (1999) con- 
ducted self-consistent N-body simulations of satellites interact- 
ing with disk galaxies and found that massive satellites may be 
accreted without destroying the primary disk galaxy, particu- 
larly for satellites on retrograde orbits which produce signifi- 
cantly less disk heating than satellites on prograde orbits. 

>From an observational point of view, Carignan et al. (1997) 
found that seven satellites around the giant lenticular galaxy 
NGC5084 exhibit retrograde motion while only one satellite 
is prograde. Moreover, the large mass and the tilted disk of 
NGC5084 suggest that this galaxy has survived a number of 
mergers, possibly with prograde satellites. Similarly, Zaritsky 
et al. (1997b) examined the orbits of their sample of satellites 
about disk galaxies and found that a slight excess (52%) of 
satellites exhibit retrograde motions in a sample of 69 primaries 
with 1 15 satellites. 

Any statistical work on satellite galaxies has to face the prob- 
lem of the small average number of satellites that can be iden- 
tified for each primary (e.g. Zaritsky et al. 1997b; McKay et al. 
2002; Prada et al. 2003). Provided that the primary galaxies 
are selected in a consistent and homogeneous way, it is useful 
to consider all satellites in the sample as satellites of a single, 
fictitious primary. This is the assumption that we work under 
here. 

In this paper we use a sample of hosts and satellites extracted 
from the SDSS (Prada et al. 2003) and a sample of halos and 
subhalos in a cosmological A^-body simulation in order to ad- 
dress a number of questions concerning the angular momenta 
of disk galaxies relative to their satellites as follows. 

1 . What is the distribution of angular distances of satellites 
from the planes of the primary disk galaxies and how does 
this compare to the theoretically-predicted angular distri- 
bution of satellite halos in simulations of structure forma- 
tion in the context of CDM? 

2. What is the ratio between prograde and retrograde satel- 
lites (P/R) in a well-defined sample of observed satellite 



systems and how does this compare to predicted satellite 
orbits? 

3. Does the fraction of prograde satellites depend on the an- 
gular distance of the satellites from the plane of the pri- 
mary disk or with projected physical separation? 

4. How does the presence of interlopers (objects that are 
field galaxies, not dynamically bound to the primary, but 
counted as satellites due to projection) affect P/R in our 
observational sample relative to the value for "true" satel- 
lites? 

This manuscript is organized as follows. In the next section, we 
describe the numerical simulations and analysis methods that 
we use to compute theoretical predictions for satellite distri- 
butions and present the basic results from the simulation. We 
describe our data and present the observational results in § 3. 
In § 4, we discuss the effects of interlopers and observational 
biases. Finally, we present our results and draw our conclu- 
sions in § 5. The main sample from which our objects are 
extracted (see Prada et al. 2003) was compiled using B mag- 
nitudes, therefore, throughout this paper, we refer exclusively 
to B magnitudes, either apparent or absolute. 

2. THEORETICAL PREDICTIONS 
2.1. Numerical Simulation 

To compare the observed satellite sample to the theoretically- 
predicted distributions of satellites in the context of CDM, 
we analyzed a simulation of structure formation in the so- 
called concordance cosmological constant plus CDM cosmol- 
ogy (ACDM), with n M = 0.3, fl A = 0.7, h = 0.7, n B h 2 = 0.22, 
« s = 1, and a power spectrum normalization of erg = 0.9. The 
simulation followed the evolution of 512 3 particles in a compu- 
tational box of 80/r'Mpc ~ 1 14.3 Mpc on a side, implying a 
particle mass of m p ~ 3.16 x 1O 8 /T 1 M . 

The simulation was performed using the Adaptive Refine- 
ment Tree A^-body simulation code (ART Kravtsov et al. 1997; 
Kravtsov 1999). The ART code employs a particle-mesh ap- 
proach and achieves high force resolution by refining the grid 
in all high-density regions with a recursive refinement algo- 
rithm. This creates an hierarchy of refinement meshes of dif- 
ferent resolution covering regions of interest. The mesh cells 
were refined up to a maximum of eight levels, corresponding to 
a minimum cell size of h pe . d ^ ~ 1 .22/r'kpc. 

We identified halos and subhalos using a variant of the Bound 
Density Maxima algorithm (Klypin et al. 1999). We began by 
computing the local density at each particle using a smoothing 
kernel of 24 particles and identified local maxima in the density 
field. We proceeded through all density peaks, beginning with 
the highest density peak and moving toward lower density, and 
marked each peak as a potential halo center. We surrounded 
each such peak with a sphere of radius rg nc | = 20/T 1 kpc and 
excluded all particles within this sphere from further considera- 
tion as a potential halo center. The parameter rg nc | is determined 
by the size of the smallest halos that we aim to identify robustly. 
After identifying potential halo centers, we iteratively removed 
unbound particles. 

We used the remaining bound particles to compute the cir- 
cular velocity profile of the halo Vd rc = y/GM{< r)/r, and the 
maximum circular velocity V max . We defined the virial masses 
and radii of halos corresponding to a fixed, spherical overden- 
sity of A = 180 times the mean density of the universe pm, so 
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that Miso = 47r(Ap M )^i8o/3 1 ■ For satellite halos or subha- 
los that are self-bound substructures within the virial radius of 
a larger parent or host halo, the outer boundary is somewhat 
ambiguous. We adopted a truncation radius r u at which the 
slope of the density profile became greater than a critical value 
of dlnp/dlnr = -0.5. This criterion was based on the fact that 
we do not expect density profiles of CDM halos to be shallower 
than this and, empirically, this definition is approximately equal 
to the radius at which the background density of the host halo 
particles is equal to the density of the particles bound to the 
subhalo. In what follows, we choose to quantify the size of ha- 
los and subhalos according to the peak circular velocity V ma x, 
because for subhalos this quantity is measured more robustly 
and this quantity is not subject to the same ambiguity as any 
particular mass definition. 

2.2. Analysis Methods 

One of our goals is to compare the observed statistics of satel- 
lite galaxies to theoretical predictions using the catalogs of ob- 
served satellites and primaries described in § 3. In order to do 
this, we construct catalogs of galaxies with absolute B mag- 
nitudes from the halos and subhalos of the simulation. There 
are several empirical methods to connect galaxies to halos that 
have been explored in the recent literature (e.g. Seljak 2000; 
Berlind & Weinberg 2002; Yang et al. 2003; van den Bosch 
et al. 2003a,b; Zheng et al. 2004; van den Bosch et al. 2005). 
For simplicity, we assign magnitudes to halos using the pre- 
scription advocated by Kravtsov et al. (2004a) in their study of 
galaxy and halo clustering statistics. We map magnitudes onto 
halos by requiring the number density of halos with maximum 
circular velocities greater than some V max to be equal to the ob- 
served number density of galaxies with absolute B magnitudes 
greater than some corresponding M B (V max ). To compute the 
number densities of galaxies of a given magnitude, we used the 
Schechter (1976) function fits to the SDSS g-band luminosity 
function presented in Blanton et al. (2003). To convert from g 
to B, we used the conversions compiled by M. Blanton at URL 
http : / / cosmo . nyu . edu/blanton/kcorrect/. The 
resulting correspondence between halo V max and M B is shown 
in Figure 1, where we also show the number density of halos as 
a function of V max . Kravtsov et al. (2004a) showed that a sim- 
ilar mapping of galaxies onto halos using r-band luminosities 
reproduces the clustering statistics of galaxies as a function of 
luminosity as measured by the SDSS. 

For each primary halo in the catalogs, we computed refer- 
ence vectors that allow us to orient satellites relative to halos in 
a meaningful way. To compute these vectors, we used all pri- 
mary halo particles within R < 0.3Rm) in order to mitigate the 
influence of large substructures at large halo-centric distances 
on these vectors and to focus on the properties of the inner halo 
which contains the material that collapsed earliest and where 
any disk galaxy would reside. First, we computed the principal 
axes of the triaxial primary halos using an iterative algorithm as 
described in detail in Zentner et al. (2005b, see also Dubinski 

6 Carlberg 1991; Kazantzidis et al. 2004a). We also computed 

7 Note that other conventions that appear in the literature define halo mass 
with respect to a virial overdensity motivated by spherical top-hat collapse (e.g. 
Lacey & Cole 1993; Eke et al. 1998), which gives an overdensity A ~ 178 in 
models with Qm = 1 and an overdensity that varies with redshift in models with 
Cl M < 1. In the ACDM model that we study, A ~ 340 at z = and A -> 178 
at high redshift (Bryan & Norman 1998, for a useful fitting function for A(z) 
see). The definition A = 180 is convenient due to the universality of the mass 
function according to this definition (see Jenkins et al. 2001, for details) 
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FIG. 1. — Assigning magnitudes to halos in the JV-body simulation. The 
horizontal axis is the maximum circular velocity of halos and subhalos V max . 
The decreasing function of Vmax is the number density of halos as a function 
of maximum circular velocity, rc(V ma x). This line should be read against the 
vertical axis on the right side of the plot. The increasing function of Vmax is 
the absolute B magnitude Mb, of galaxies with the same number density as 
the number density of halos at this value of V max (decreasing line). This line 
should be read against the vertical axis on the left side of the plot and it repre- 
sents the mapping of halo size (Vmax) to absolute magnitude (Mb) that we use 
to compare our observational results to halos and subhalos in dissipationless 
simulations. 



the net angular momentum vector of all primary halo particles 
within R < 0.3Rm)- In the following, we compare various hy- 
potheses for the orientations of disks in their primary halos to 
the data regarding the relative positions and velocities of satel- 
lite galaxies in our sample. We computed all statistics by sum- 
ming over three orthogonal projections through the simulation 
volume and we computed statistical errors by taking each pro- 
jection to be statistically independent. 

2.3. Theoretical Predictions 

After assigning magnitudes to halos and determining refer- 
ence vectors for halos, we constructed catalogs of primaries 
and satellites by observing the simulation from three orthog- 
onal projection directions using the same criteria as we em- 
ployed for our observational sample (see § 3.1 below). Pri- 
maries were restricted to only those halos that are assigned 
magnitudes -20.5 < Mb < -19.5. Primaries were required to 
satisfy isolation criteria such that any halos with projected sep- 
arations r p < 500/r'kpc and line-of-sight velocity differences 
AV < 1000 km s" 1 with respect to the primary are dimmer by 
at least AM B > 2. Satellites were required to have AMb > 2, 
r p < 350/z _1 kpc, and AV < 500 km s" 1 . 

The first results from the simulations concern the angular 
positions of the observed satellites. We have examined three 
idealized hypotheses regarding the orientation of disks in dark 
matter halos: (1) that the disk rotation axis is aligned with the 
angular momentum of the halo; (2) that the disk is aligned with 
the major axis of the host halo; (3) that the disk is aligned with 
the minor axis of the halo (which is generally in close alignment 
with the angular momentum axis). The first hypothesis is based 
recent attempts to reconcile the anisotropic of distribution of 
satellites in the Local Group (e.g. Lynden-Bell 1982; Hartwick 
2000; Willman et al. 2004; Kroupa et al. 2005) with the pre- 
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dictions of numerical simulations (e.g. Zentner et al. 2005a; 
Libeskind et al. 2005). Alternatively, the simplest picture of 
disk galaxy formation predicts that halo and disk angular mo- 
menta should be nearly aligned (Fall & Efstathiou 1980, though 
cosmological hydrodynamic simulations show a misalignment 
of baryon and dark matter angular momenta that is typically 
- 10° -20°, see van den Bosch et al. 2002; Chen et al. 2003). 
We have defined the angle cf) as the angle between the long 
axis of the primary disk and the position of the satellite in two- 
dimensional projection. 

The cumulative angular distributions are shown in Figure 2. 
In the case where the disk is closely aligned with the major 
axis of the primary halo (solid lines), the anisotropy is appar- 
ent as a deficit of satellites at small and intermediate angles in 
Figure 2. In this scenario, satellites are located preferentially 
near the short axes of their disk primaries in a sense similar 
to that observed by Holmberg (1969), Zaritsky et al. (1997a), 
and Sales & Lambas (2004). A Kolmogorov-Smirnov (KS) 
test reveals that the probabilities for the simulated samples with 
Mb < — 16.5 and Mg < -15.5 to be drawn from an isotropic un- 
derlying distribution are P K s — 9 x 10" 5 and P K s — 5 x 10" 5 
respectively. In Figure 2, we also show the resulting distribu- 
tions under the assumption that disk rotation is aligned with the 
angular momentum axis of the inner halo of the primary. Un- 
der this assumption, their appears to be a very slight excess of 
satellites at small and intermediate angles, but the distributions 
are consistent with isotropy. The distributions with the disk an- 
gular momenta aligned with halo minor axes are quite similar 
to these and we have omitted them from Figure 2 in the interest 
of clarity. 

The next results concern the fraction of satellites that exhibit 
prograde or retrograde motion with respect to the disk of the 
primary galaxy. In the case of the alignment of disk rotation 
with the major or minor axes of the host halo, there is no obvi- 
ous way to assign a sense to the rotation to the disk. For this rea- 
son, we have computed prograde and retrograde statistics only 
for the case where the direction of the disk rotation is assumed 
to be aligned with the net angular momentum vector of the dark 
matter of the primary halo within R < 0.3Rm)- The results are 
summarized in Table 1. The error estimates listed in Table 1 
represent the standard deviation of a binomial distribution with 
the appropriate parameters for each case. By presenting the 
prograde/retrograde statistics in this way, we have assumed that 
we can treat the three projections through the simulation as in- 
dependent. The excess of prograde satellites is a robust predic- 
tion in this scenario. Additionally, notice the trend for larger 
satellites to have a greater probability to exhibit prograde rela- 
tive motion. For example, the satellites in the largest bin show 
a significant trend toward prograde motion as do the satellites 
with -17.5 > M B > -18.0, where there are 247 satellites with 
141 progrades, yielding f plog =0.57 ±0.03. Alternatively, in the 
magnitude range -15.5 > M B > -16.0, there are 294 satellites, 
152 of which are prograde, yielding a prograde fraction in this 
bin of luminosity of f pmg = 0.52 ± 0.03. 

The excess of prograde satellites and the trend toward higher 
prograde fractions for larger satellites are not altogether sur- 
prising. CDM halos tend to be connected through a network of 
filamentary structures established by the statistics of the den- 
sity field and the tidal field during collapse in overdense re- 
gions (e.g. Klypin & Shandarin 1983; Bond & Myers 1996; 
Bond et al. 1996; Colberg et al. 2004). The matter distribution 
within the filaments is strongly concentrated toward the axis of 
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FIG. 2. — The cumulative angular distribution of satellites relative to the 
long axes of their primaries in simulations. The three panels correspond to 
three different magnitude limits, Mb < -17.5, Mb < —16.5, and Mb < -15.5 
from top to bottom. Solid lines correspond to the assumption that the disk 
rotation axis is aligned with the major axis of the primary halo. Dashed lines 
correspond to the assumption that the disk rotation axis is aligned with the 
angular momentum of the inner halo of the primary. The straight, dotted lines 
correspond to an isotropic distribution in <j>. In the lower right of each panel, 
we give the KS probability for the satellites in the samples with disk rotation 
along the major axis of the halo to be drawn from an isotropic distribution. In 
all cases, the simulated samples with disk rotation aligned with halo angular 
momentum are consistent with isotropy. 



the filament (Colberg et al. 2004) and halos preferentially tend 
to form in these dense regions and then merge along the dom- 
inant filamentary directions. This leads to a strong correlation 
between directions along which halos merge during the mass 
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Table 1 

Prograde and Retrograde Satellites In Simulations 
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The prograde and retrograde fractions of satellites in the simulation 
volume. Column description: (1) The minimum value of V max for satel- 
lites in the sample in kms -1 ; (2) The corresponding maximum magni- 
tude determined by matching halo number densities to galaxy number 
densities; (3) The total number of satellites in the sample, summed 
over the three, orthogonal projections through the simulation volume; 
(4) The number of prograde satellites that have peculiar velocities that 
are in the same direction as the rotation of the inner halo; (5) The num- 
ber of retrograde satellites that have peculiar velocities that are in the 
opposite direction of the rotation of the inner dark matter halo; (6) The 
fraction of satellites with prograde velocities. 



accretion history of the primary halo (Knebe et al. 2004; Zent- 
ner et al. 2005b). The subhalos that merge early and establish 
the angular momentum of the inner regions of the primary halo 
(e.g. Vitvitska et al. 2002) flow along directions that are similar 
to the directions of infalling subhalos at late times, which are 
then observed as the satellites. Larger halos tend to be more 
strongly biased toward formation in the overdense regions (e.g. 
Bond & Myers 1996; Bond et al. 1996), so they are more faith- 
ful tracers of the flow along the filaments. The higher prograde 
fractions for larger satellites are a reflection of this. 

The excess of prograde satellites is also shown in Figure 3, 
which is a histogram of the peculiar velocities of the simulated 
satellite halos with respect to their primaries. In this plot, pos- 
itive peculiar velocities are defined to be in the same sense as 
the rotation of the disk. Again, the modest excess of prograde 
satellites, particularly with relatively small peculiar velocities, 
is apparent. 

The simulated catalog provides a tool to study the effect of 
interlopers on the prograde/retrograde fractions because in this 
case we know both the "observed" sample and the true three- 
dimensional separation. We have taken as a working definition 
to separate true satellites from interlopers a three-dimensional 
distance of r true = 1 /r'Mpc. Satellites that are a distance greater 
than r tme from their primaries are considered interlopers. The 
virial radii of primaries in our sample range from roughly 
175/r'kpc to 400/r'kpc, so we generally expect objects at dis- 
tances greater than ~ 1 /z _1 Mpc to be unrelated to the primary. 
In Table 2, we show the prograde/retrograde fractions for the 
true satellites separated in three-dimensions by a distance less 
than r true , with no contamination by interlopers. We find that 
the net effect of interlopers at large line-of-sight separations is 
to dilute the prograde/retrograde fractions by ~ 1 - 7%, depend- 
ing upon the details of the sample selection. 
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FIG. 3. — The peculiar velocities of the simulated satellite samples relative 
to the angular momenta of the central regions of the primary host halos. The 
histograms show the fraction of satellites in bins of peculiar velocity with ret- 
rograde satellites on the left (negative peculiar velocity relative to the disk) 
and prograde satellites on the right. The vertical, dashed line shows the sepa- 
ration between prograde and retrograde satellites at zero peculiar velocity. We 
show histograms for three subsamples labeled by both their minimum Vmax 
and their maximum Mb . The solid line IS Vmax >76kms-' (M B < -15.5), the 
dotted line is for Vmax > 82kms~' (Mb < -16.5), and the dashed line is for 
Vmax > 92kms-' (M B < -17.5). 

Table 2 

Prograde and Retrograde Satellites In Simulations 
with Interlopers Removed 





Afprog 


Nretr 


/prog 


/false 


-18.0 


42 


16 


0.72 ±0.06 


0.29 


-17.5 


145 


90 


0.62 ±0.03 


0.29 


-17.0 


278 


196 


0.59 ±0.02 


0.27 


-16.5 


385 


291 


0.57 ±0.02 


0.26 


-16.0 


499 


389 


0.56 ±0.02 


0.27 


-15.5 


608 


496 


0.55 ±0.02 


0.27 



The prograde and retrograde fractions of satellites with true, three- 
dimensional distances smaller than r true = l/r'Mpc from their asso- 
ciated primaries in the simulation volume. Column description: (1) 
The corresponding maximum magnitude determined by matching halo 
number densities to galaxy number densities; (2) The number of pro- 
grade satellites that have peculiar velocities that are in the same di- 
rection as the rotation of the inner halo; (3) The number of retrograde 
satellites that have peculiar velocities that are in the opposite direction 
of the rotation of the inner dark matter halo; (4) The fraction of pro- 
grade satellites; (5) The fraction of all satellites in Table 1 that were 
identified as "false" satellites at each magnitude threshold, defined by 
the fact that they were included in the sample due to projection effects 
but are a distances greater than r trac = 1 /r'Mpc from their primaries. 



3. OBSERVATIONAL RESULTS 



3.1. The Observed Sample 
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Our observed sample of systems is extracted from "Sample 
2" of Prada et al. (2003), which was constructed from SDSS 
data. Sample 2 has a maximum depth of 60000kms _1 and 
a limiting absolute magnitude in B of Mb = -18. The isola- 
tion criteria of the primaries in the sample are that any satel- 
lites within a projected radial separation of r p < 500/r'kpc 
and with a line-of-sight velocity difference AV < lOOOkms" 1 
must be dimmer than the primary by at least AMb > 2. Satel- 
lites must satisfy AMb > 2 with respect to their associated 
primary, they must have a projected distance to the primary 
of r p < 350/i _1 kpc, and a line-of-sight velocity difference of 
AV < 500 km s" 1 . 

Sample 2 contains more than 700 primaries and more than 
1000 satellites. Using Sample 2 as a starting point, we selected 
our primaries according to the following additional criteria. 

1 . We enforced a limiting depth of 1 1000 kms" 1 . 

2. We chose primaries only in the magnitude range -20.5 < 
M B < -19.5 

3. We required all primaries to have a well-defined morpho- 
logical type as a spiral. 

The limit of 11000 kms" 1 was imposed not to lose too many 
satellites, as the SDSS is complete down to Mb = -17.5 at 
11 000 kms" 1 , thus the satellites (two magnitudes fainter than 
primaries) are those brighter than Mb = -15.5 only. The narrow 
magnitude range for the primaries was chosen to avoid includ- 
ing primaries with very different masses, given the known de- 
pendence of luminosity on the mass of the primary (e.g. Prada 
et al. 2003). Further, satellites with projected distances from 
their primaries smaller than 20kpc were rejected, due to po- 
tential galaxy misclassification (e.g. HII region contamination, 
bright stars etc.). After applying these constraints, we were left 
with 144 spiral primaries and 193 satellites. We obtained the 
direction of rotation of a limited number of primaries, so we 
analyzed the prograde and retrograde motion using a subsam- 
ple of 43 primary galaxies with 76 satellites (see § 3.3). For the 
angular distribution analysis, we used the entire sample of 193 
satellites (see § 3.2). 

The observational data listed in Table 3 come from the SDSS 
database and from our own observations with the exception 
of NGC2841, which had a previously -published rotation curve 
that we took from Afanasiev & Silchenko (1999). Data such 
as recessional velocities, absolute magnitudes and positions for 
primaries and satellites were taken from the SDSS database, 
while the Position Angles of the primaries were taken from 
the Lyon-Meudon Extragalactic Database (LED A) 8 . Low- 
resolution spectroscopy of the Ha line easily determines the 
spin direction of our objects. Therefore, through an observa- 
tional program at the telescopes of the Isaac Newton Group (La 
Palma, Canary Islands), we obtained rotational data for 42 of 
our primaries. 

Observations were performed using the spectrographs ISIS, 
of the 4.2 m William Herschel Telescope (WHT), and IDS, of 
the 2.5 m Isaac Newton Telescope (INT) of the Isaac Newton 
Group, on La Palma, Spain. The setup was simple and no spe- 
cial observing conditions were needed, so most of the data were 
taken in service mode. The most important aspect of the obser- 
vation was to determine reliably the orientation of the charge- 
coupled device (CCD), to be able to place the receding end of 

8 URL http://leda.univ-lyonl.fr 



0.8 



0.6 

V 



0.4 



0.2 





1 1 1 1 1 1 1 1 1 1 1 
Whole sample, 193 satellites 


1 ff 




Vol-lim. (M<-17.5), 74 sat. 




- 


Vol-lim. (M<-16.5), 64 sat. 




- 


rX 

s'f/ 

f 

S 










- 


yd-''" 










Jl / 






/■*•' 




- 


&/ 


- 




f/'t 
















i i 1 i i i 1 i i i 1 i i 





20 40 60 80 

[°] 

FIG. 4. — The cumulative distribution of the angular positions of satellites 
relative to their primary disk galaxies. The samples shown in this plot are 
the entire sample with 193 satellites and 144 primaries (solid line), a volume- 
limited subsample at a depth of 11000 km s with 74 satellites and 62 pri- 
maries (dashed line) and a volume-limited subsample at a depth of 7000 km s~ 1 
with 64 satellites and 51 primaries (dot-dashed line). The thin, dotted line in 
the background represents an isotropic cumulative distribution function. 



the galaxies correctly on the sky. Two exposures at Sky PA 0°, 
shifted by some 10", and two more at Sky PA 90° with the same 
offset were sufficient to mark North and East on the images. 

We used the R316R or R600R gratings on the Marconi II 
chip with ISIS, centered at 6650, and R900V, R600V on chip 
EEV13 or R1200R on chip EEV10 with IDS, also centered at 
6650. The slit width ranged from 1" to 1.5", depending on the 
seeing conditions. Exposure times were typically 300 s with 
WHT and 900 s with INT. During the observations, the seeing 
varied from 0.8" up to 2.5". The detection of the shift on the 
Ha line was performed by analyzing the images with the IRAF 
package, comparing the line centroid shift with that of a sky 
line to correct for spurious rotations of the CCD. The observing 
runs took place over the period from March to September of 
2003. 

3.2. The Angular Distribution of Satellites 

In order to study the distribution of the angular positions of 
the observed satellites relative to their disk primaries we were 
able to use our largest sample of objects because this requires 
only the coordinates of all objects and the Position Angles of 
the primaries, which we extracted from the LEDA database. 
We included in this analysis the 144 disk primaries in the mag- 
nitude bin -20.5 < M B < -19.5, and the 193 satellites selected 
according to the description in § 3.1. We reduced the angular 
distance of the satellites from the disks to a single quadrant, 
0e[O°,9O°]. 

Figure 4 shows the cumulative frequency for the entire satel- 
lite sample (solid line), compared to that of an isotropic distri- 
bution (dotted line). In addition, we computed the angular dis- 
tributions of satellites in two volume-limited subsamples. The 
first subsample has a depth of 1 1000 kms" 1 and limiting magni- 
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Table 3 

The Primaries and Satellites of the sample. 



1 

Name 


Bmag 2 


RA(2000) ' 


DEC(2000) ' 


— — a — 

PA 4 


V(kms _1 ) s 


1 T. 

Ay 6 


dist (kpc) ' 


Orbit s 


pgc001112 


-20 


.3 


00 


16 


55 


-00 05 18 


170 


10042.9 


- 


- 


- 


a 


-18. 


.3 


00 


17 


07 


-00 08 30 


137 


10098.7 


55.8 


175.1 


r 


pgc001841 


-20 


.1 


00 


30 


07 


-11 06 49 


169 


3649.6 








a 


-15. 


.9 


00 


31 


28 


-10 40 33 


37 


3645.2 


-4.4 


486.0 


P 


mcg-2-3-61 


-20. 


,2 


01 


00 


04 


-11 04 56 


123 


5519.4 








a 


-17. 


,9 


01 


00 


05 


-11 02 32 


6 


5420.0 


-99.4 


52.2 


r 


ngc341 


-20. 


.4 


01 


00 


45 


-09 11 08 


55 


4658.0 








a 


-16. 


.9 


01 


01 


43 


-09 19 00 


119 


4639.3 


-18.7 


294.5 


P 


ugcl962 


-20. 


.0 


02 


28 


54 


00 22 13 


105 


6152.9 








b 


-17. 


.4 


02 


29 


33 


00 22 23 


91 


6549.3 


396.4 


240.4 


P 


c 


-17. 


.0 


02 


29 


39 


00 07 23 


37 


6155.4 


2.5 


455.2 


P 



Host names and satellite labels 
2 Absolute magnitude in B 
3 J2()()0 coordinates of the objects 

4 Sky PA of the hosts and, for the satellites, the Position Angle of their location with respect to the host 

5 Recessional velocity (h = 0.7) 

6 The difference in recessional velocity 
'Projected distance between satellite and host 
8 "p" stands for prograde and "r" for retrograde. 

Note. — This table is published in its entirety in the electronic edition of the Astrophysical Journal. 
A portion is shown here for guidance regarding its form and content. 



tudeMfi =-17.5, and the second has a depth of 7000 kms" 1 and 
limiting magnitude M B = -16.5. The statistics for the volume- 
limited subsamples are poorer (74 satellites at 1 1000 kms" 1 and 
64 at 7000 kms -1 ), but the behavior is qualitatively similar to 
that of the entire sample of 193 satellites. 

At first glance, this result seems most easily compatible with 
the simulated distribution with disk rotation nearly aligned with 
the angular momentum of the inner halo (Figure 2). We have 
performed a comparison of the distributions, using a KS test 
to test several hypotheses. The full observational sample as 
well as the volume-limited subsamples are all consistent with 
an isotropic distribution about their primaries. Contrary to the 
results of Zaritsky et al. (1997a) and Sales & Lambas (2004), 
this sample does not show significant evidence of anisotropy. 
In addition, the observed sample is consisted with all three of 
the hypotheses used to construct distributions from the simula- 
tion. This includes the case with the disk rotation axis aligned 
with the major axis of the host halo, which yields a probability 
Pks — 0.17 of being drawn from the same underlying distribu- 
tion. The relatively small size of the current observational sam- 
ple makes it difficult to reject any of our idealized disk align- 
ment hypotheses. 



3.3. Observed Prograde Satellite Fractions 

To state if a satellite is prograde or not, we need to know the 
difference in recessional velocity between the satellite and the 
primary (Av = V sa , - V pr j m ), the Position Angle of the satellite 
with respect to the primary (PAS), and the sky Position Angle 
of the redshifted end of the primary disk (PAred). We define 



Table 4 

Prograde/Retrograde Statistics of the Entire 
Observational Sample 



Angle from disk 


Mats 


"prog 


Nretr 


/prog 


0° to +30° 


24 


15 


9 


0.63±0.10 


0° to +45° 


40 


25 


15 


0.63 ±0.08 


0° to +75° 


58 


38 


20 


0.66 ±0.06 


0° to +90° 


76 


46 


30 


0.61 ±0.05 



Statistics of prograde/retrograde satellites for the entire sample of 
satellites. The columns are (1) the angular displacement from the pri- 
mary disk plane, summed over the four quadrants; (2) the number of 
satellites for each angular selection; (3) the number of prograde satel- 
lites found for each angular selection; (4) the number of retrograde 
satellites found for each angular selection; (5) the frequency of pro- 
grade satellites over the total with an error estimate given by the stan- 
dard deviation of the appropriate binomial distribution. 

the quantity S by 

J | PAred- PAS |< 90° \ 
i -+l tt <jj pAred _ pAS j> 270° J 

S = -l otherwise, 

A satellite is prograde if S ■ Av > 0, and retrograde S ■ Av < 
0. Table 3 gives the sense of the line-of-sight velocities of the 
satellites in our sample with respect to their disk primaries. 

The resulting prograde satellite fractions are summarized in 
Table 4 and Table 5. Unfortunately, the small size of the ob- 
servational sample does not allow us to study any magnitude 
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Table 5 

Prograde/Retrograde Statistics for the Entire Sample 
to a Depth ofV< 7000km s _1 (upper panel) and for a 
Volume-limited Sample to a Depth ofV< 7000km s -1 

(LOWER PANEL) 



Angle from disk 






AW 


/prog 


0° to +30° 


11 


8 


3 


0.72±0.13 


0° to +45° 


20 


14 


6 


0.70±0.10 


0° to +75° 


34 


22 


12 


0.64 ±0.08 


0° to +90° 


41 


26 


15 


0.63 ±0.08 


Angle from disk 


AU 


"prog 


Aretr 


/prog 


0° to +30° 


9 


6 


3 


0.66±0.16 


0° to +45° 


13 


8 


5 


0.61 ±0.13 


0° to +75° 


23 


13 


10 


0.57±0.10 


0° to +90° 


28 


16 


12 


0.57 ±0.09 



Prograde/Retrograde statistics of satellite galaxies for an unlimited 
subsample at a depth of 7000km s -1 and for a corresponding volume- 
limited subsample at the same depth. The depth of 7000 km s -1 (Mb = 
-16.5) has been chosen in an attempt to keep a significant number of 
object in each subsample. The columns are (1) the angular displace- 
ment from the primary disk plane, summed over the four quadrants; 
(2) the number of satellites for each angular selection; (3) the number 
of prograde satellites found for each angular selection; (4) the num- 
ber of retrograde satellites found for each angular selection; (5) the 
frequency of prograde satellites over the total. 
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FIG. 5. — The distribution of peculiar velocities of the entire unlimited 
sample (solid line, empty), and the volume-limited sample at a depth of 
7000 kms -1 (dotted line, dashed). Positive peculiar velocities correspond to 
prograde satellites and negative peculiar velocities correspond to retrograde 
satellites. The zero peculiar velocity is marked by a vertical dashed line. 



dependence of the prograde satellite fraction at a statistically- 
significant level as we did for the simulated satellites in § 2.3. 
Alternatively, we present the prograde/retrograde fractions as 
a function of angular displacement from the major axis of the 
primary galaxy disk. Table 4 shows the resulting distribution of 
prograde and retrograde satellites for the entire sample at full 
depth (1 1000 kms -1 ). Table 5 shows two subsamples with lim- 
iting depths of 7000 kms" 1 , one containing all objects to this 
depth and one for a volume-limited subsample with an absolute 
magnitude limit of Mb = -16.5. 

It is evident from Table 4 that the largest samples of satel- 
lites exhibit an excess of prograde peculiar velocities that is 
significant at the ~ 2 -2.5c level. In addition, the data show 
no evidence for any significant trend in the prograde fraction as 
a function of the angular position with respect to the disk. At all 
angles with respect to the disk, the prograde fraction is ~ 60%, 
though the statistical significance of the prograde excess in the 
smallest angular cuts is marginal. Table 5 shows a qualitatively 
similar result, with f pmg ~ 60% in both of these subsamples, 
though again, the statistical significance is, at best, marginal in 
these small samples. In addition to these results we also found 
no significant evidence for a variation in the prograde fraction 
as a function of projected distance to the primary. 

Figure 5 shows a histogram of the peculiar velocities of 
two of our samples, the entire sample with all objects and the 
volume-limited subsample at a depth of 7000kms _1 (M B < 
-16.5). Figure 5 shows a peak in the slowly-rotating progrades 
in both of these samples. This seems to be in broad agreement 
with the results of the simulated sample of subhalos shown in 
Figure 3. A more detailed analysis would require a larger ob- 
servational sample. 

4. THE EFFECT OF OBSERVATIONAL BIASES 



We performed a series of Monte Carlo simulations of mock 
observational samples in order to quantify the effects of obser- 
vational biases and interlopers. We produced lists of objects, 
both satellites and interlopers, distributed in magnitude accord- 
ing to the luminosity function of each class of object. This anal- 
ysis complements the previous results on the net effect of in- 
terloper contamination in the cosmological /V-body simulation 
performed in § 2.3), because it allows us to study a statistically- 
large number of independent, mock observations for a variety 
of assumptions about the true fraction of prograde satellites. 

The first set of Monte Carlo simulations we performed were 
aimed at assessing the effect of interlopers (which are sepa- 
rated from their primaries by large distances along the lines 
of sight, but are associated due to projection) on the measured 
prograde/retrograde fractions of satellites relative to the pro- 
grade/retrograde fractions of "true" satellite galaxies that have 
small three-dimensional separations from their primaries. The 
number of objects at each bin of magnitude was determined 
using Schechter (1976) luminosity functions with different pa- 
rameters according to each class of object. Specifically, we em- 
ployed the magnitude form of the Schechter function, 

$(M) = 0.4$»log(10) io°- 4 ( M *- M X« +1 ) e xp(-10 ' 4(M -- M) ), (1) 

which we then used to determine the number of objects by in- 
tegrating over magnitude and volume. 

For the interlopers, we used the parameters given by Blan- 
ton et al. (2003) for SDSS field galaxies, converted to B band 
(see section 2.2), namely = 0.0074 Mpc" 3 , M* = -21.31, 
a = -0.89. In this case, the volume we considered was the 
same truncated cone as in our real observations, but with an ad- 
ditional constraint that the distance r from the primary satisfies 
r > 1 /i _1 Mpc. This constraint serves as our working definition 
of "interloper" as opposed to "true" satellite. 
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For the satellites, we used the parameters of the Local Group 
(LG), assuming that it represents a typical system of satellite 
galaxies (Mateo 1998; van den Bergh 2000). The Schechter 
function parameters for the LG satellites were M* = -19.5 and 
a = -1.14. However, the characteristic number density <£>™ f 
for the LG is not well defined as it depends sensitively on the 
volume considered. A way to get a first approximation of its 
value is to assume a given interloper contamination, which can 
be taken from our Table 1 in § 2.3, and is in agreement with 
the fraction of interlopers found by Prada et al. (2003) from 
their satellite dynamics analysis. We then assume an interloper 
contamination of 25% and calculate $ f such that we recover 
25% of interloper contamination in the simulation. We used 
$'«*' ^ 0.07 Mpc~ 3 . The volume considered for the satellites 
was the intersection of the truncated cone of our real observa- 
tions with the sphere contained within r < 1 /r'Mpc about each 
primary. 

We then observed the objects in our simulation in the same 
way as we did for the SDSS data. For simplicity, in our first 
suite of Monte Carlo mock observations, we kept a fixed re- 
cessional velocity of 7000kms _1 for all the objects in these 
mock catalogs. We assumed that interlopers always have a 
50% chance of being prograde or retrograde, while the prograde 
fraction of satellites and the total number of primaries are both 
input parameters that are used to construct specific mock obser- 
vations. The number of satellites for each galactic system (or 
for each primary) is then calculated using the Schechter func- 
tion with the <f>'™ f as above. 

We run 500 Monte Carlo trials, with an input prograde frac- 
tion of / prog = 60%; the input number of primary galaxies was 
selected such that the total number of combined false and true 
satellites in the mock sample was about 76 as in our real ob- 
servations. The interloper contamination in our mock sam- 
ples reduced the observed prograde fraction to an average of 
/prog = 0.57 with a root-mean-square scatter about the mean of 
±0.05. This agrees quite well with the effect of interlopers 
found in the mock catalogs constructed from the halos in the 
cosmological A^-body simulation (see Tables 1 and 2). 

The second set of Monte Carlo simulations was aimed at 
determining how the prograde/retrograde fractions are affected 
by having systems at different recessional velocities mixed to- 
gether, as we have in our observed sample. When several 
galactic systems at different recessional velocities are mixed 
together, different levels of interloper contamination are super- 
posed. Also, the apparent magnitude limit of the observations 
results in a selection of satellite luminosities that depends on 
the recessional velocities of the systems. The aim is to simu- 
late these two observational effects. We used the same setup as 
for the Monte Carlo simulations with fixed recessional veloci- 
ties that we discussed in the previous paragraph. We built up 
three samples of systems, again with the necessary input num- 
ber of primaries so as to "observe" a total of 76 false and true 
satellites combined, at recessional velocities of 2800,5600 and 
8800 kms" 1 . We then observed all of the samples together as a 
single, mock sample. 

We ran a series of 500 such random trials with an input satel- 
lite prograde fraction of / prog = 60% as above. The mean ob- 
served percentage of prograde satellites in the mock samples 
was /p r0 g ~ 58%, with a root-mean-square scatter of ~ 5% 
about the mean. This implies that the systematic dilution of the 
prograde/retrograde fraction due to interlopers is fairly small 
for comparable samples; however, the scatter about the true 




True (input) prograde fraction (%) 

FIG. 6. — Observed reduction of the prograde fraction caused by interloper 
contamination in a sample containing systems at a wide range of recessional 
velocities. The horizontal axis represents the true prograde fraction for all true 
satellites and the vertical axis represents the observed prograde fraction taking 
interloper contamination into account. The points represent the mean observed 
prograde fraction derived from our Monte Carlo simulations of the net effect 
of interloper contamination in the case of mixing systems with different reces- 
sional velocities. The error bars represent the root-mean-square scatter about 
the mean. The dotted line represents perfect retrieval of the prograde fraction. 



fractions due to the presence of interlopers is not insignificant 
for samples of this size. We show in Figure 6 the mean reduc- 
tion of the "observed" prograde fraction as a function of the true 
input prograde fraction as well as the la scatter in the observed 
values about the mean. 



5. CONCLUSIONS 

We study the distribution of satellite galaxies about disk 
primaries. The study has both observational and theoretical 
components. We focus our attention on the fraction of satel- 
lite galaxies that exhibit prograde motion with respect to their 
disk primaries and on the angular distribution of satellites rel- 
ative to their primaries. On the observational side, we study a 
carefully-selected observational sample of primary disk galax- 
ies and satellite galaxies extracted from the Sloan Digital Sky 
Survey database. In addition, we analyze the structure and 
distribution of halos and subhalos in an N-body simulation of 
structure formation in the, now standard, concordance ACDM 
model. 

We find that the fraction of prograde satellites in our entire 
flux-limited sample is / pr og = 0.61 ±0.05. The observational 
samples show no evidence for a strong dependence of /p r0 g on 
magnitude, angular distance from the satellite to the major axis 
of the primary disk, or projected distance from the primary. 
However, the observed samples are, as yet, small and as a result, 
the statistics are poor. In addition, we found that the observed 
satellite galaxies are consistent with being distributed isotropi- 
cally about their primary disk galaxies. 
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From the cosmological A^-body simulation, we constructed 
catalogs of primaries and satellites using an algorithm that maps 
galaxies onto halos and subhalos based on matching number 
density (e.g. Kravtsov et al. 2004b). We built catalogs based 
on three idealized hypotheses for the orientations of disk galax- 
ies within their surrounding dark matter halos, with the angular 
momentum of the disk aligned with either (1) the major axis of 
the inner host halo, (2) the minor axis of the host halo, or (3) the 
angular momentum of the particles in the inner host halo. We 
then presented the general predictions of these scenarios. The 
first scenario predicts a projected distribution of satellite galax- 
ies that is strongly anisotropic, with galaxies located near the 
poles of their disk primaries. For example, the Kolmogorov- 
Smirnov probability that the entire simulated distribution of ha- 
los that are assigned magnitudes Mb < -15.5 is selected from 
an underlying isotropic distribution is Pks — 5 x 10" 5 . How- 
ever, a KS test does not yet allow us to reject this hypothe- 
sis based on the present data, largely because the observational 
data set is fairly small (193 satellites in our largest, flux-limited 
sample). The other two hypothesized disk orientations pre- 
dict satellite distributions that are consistent with isotropy at 
all magnitudes. Under the hypothesis that the angular momenta 
of the disk primaries are aligned with the angular momenta of 
the inner regions of the dark matter halos that host them, the 
simulation predicts a prograde satellite fraction « 55-60%. 
Specifically, for the entire sample of satellites withM B < — 15.5, 
we find / prog = 0.54 ± 0.01, while for a subsample of the largest 
satellites with Mb < —17.5, we find / prog = 0.59 ± 0.03. As the 
previous sentence indicates, we find that / prog is a weak function 
of satellite size (magnitude), whereby larger (more luminous) 
satellites are more likely to exhibit prograde motions. The sim- 
ulation results are broadly consistent with the values measured 
in the observational sample. A larger observational sample will 
be needed in order to test for any magnitude dependence in the 
observed prograde satellite fraction and to make more mean- 
ingful comparisons between theoretical predictions and obser- 
vations. Larger observational samples will also require a sig- 
nificant refinement of theoretical predictions. 

In addition to these results, we also made an effort to assess 
the net effect of interlopers on the observed prograde satellite 
fractions compared to the prograde satellite fractions for "true" 
satellites. From our analysis of both the A^-body simulation and 
an independent set of Monte Carlo simulations of mock obser- 
vations, we found that the mean dilution of the prograde satel- 
lite fraction due to interloper contamination is of the order of a 
few percent. However, we found that the scatter in the observed 
prograde and retrograde fractions in samples comparable to the 
size of our observational set is sizable at nearly 5%. 

A comparable observational study including the pro- 
grade/retrograde satellite fractions and the angular distributions 
of satellites was performed by Zaritsky et al. (1997a,b) who had 
a slightly larger sample of systems. Contrary to our findings, 
Zaritsky et al. (1997b) reported more retrograde satellites than 
prograde satellites (/prog — 48%). However, note that this result 
is not statistically significant as an estimate of the statistical er- 
ror on the mean from a sample of the size of Zaritsky et al. 
(1997a) gives cr(f pmg ) ^ 5%. Zaritsky et al. (1997a) also re- 
ported strong evidence for an anisotropic satellite distribution 
as we mentioned in § 1 . We were unable to confirm this result 
with our observed satellite sample, which is consistent with an 
isotropic underlying distribution. The nature of any disagree- 
ments between our results and the previous work of Zaritsky 



et al. (1997a,b) is unclear at present; however, it could possibly 
be due to different criteria in the selection of the objects, un- 
accounted for some observational biases, and/or statistical fluc- 
tuations in the relatively small samples. Larger samples con- 
structed from forthcoming data sets should be able to address 
these issues as well. 
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